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We calculate the damping of low-lying collective modes of a trapped Bose gas in the hydrodynamic 
regime, and show that this comes solely from the shear viscosity, since the contributions from bulk 
viscosity and thermal conduction vanish. The hydrodynamic expression for the damping diverges 
due to the failure of hydrodynamics in the outer parts of the cloud, and we take this into account by 
a physically motivated cutoff procedure. Our analysis of available experimental data indicates that 
higher densities than have yet been achieved are necessary for investigating hydrodynamic modes 
above the Bose-Einstein transition temperature. 
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In recent experiments on magnetically-trapped atomic 
vapors, alkali atoms have been cooled to temper- 

atures at which they are degenerate and indeed Bose- 
Einstein condensation has been observed in them. Fre- 
quencies and damping rates of collective modes in these 
systems have been investigated, both above and below 
the Bose-Einstein transition temperature, T c Q-ra]. In 
this Letter we shall focus on properties above T c . One 
can distinguish two regimes, the hydrodynamic one, for 
which the characteristic mode frequency is small com- 
pared with the collision frequency and the wavelength 
of the mode is large compared with the atomic mean 
free path, and the opposite limit, the collisionless one, 
for which collisions are relatively unimportant. The fre- 
quencies of modes in the hydrodynamic regime have been 
calculated in Ref . , and here we calculate their damp- 
ing. 

We begin by giving a simple derivation of the basic hy- 
drodynamic equations. Our treatment is essentially that 
of Ref. || generalized to take into account the potential 
of the trap, and, since we are interested in small oscil- 
lations, we shall consider the linearized equations. The 
Euler equation for the fluid velocity v(r, t) is 
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mn (v)- 



- Vp(r, t) + mn(r, f)f, 



(1) 



where m is the mass of the atoms, n(r, t) is the particle 
density, no(r) is the equilibrium particle density, p(r,t) 
is the pressure and f is the force per unit mass due to 
the external potential Uo(r), f = -Vf/o(r)/m. In equi- 
librium, where the pressure ispo( r )> Eq. (pj) implies that 
Vpo(r) = mno(r)f. Taking the time-derivative of Eq. 
(m) and using the continuity equation, one finds 



, dp(r,t) 

at 



V ■ [n (r)v] mf . 



(2) 



We calculate the first term on the right hand side of Eq. 
(0) by using the energy conservation condition M, 



dt 



(pe) = - V • (wprv) + pv ■ f, 



(3) 



where p is the mass density and e and w are, respectively, 
the internal energy and the enthalpy of the fluid per unit 
mass. Since we assume local thermodynamic equilibrium 
and neglect contributions to the energy due to interpar- 
ticle interactions, we may use the results p — p(w — e) 
and pe = 3p/2, and find that 



dp(r,t) 
dt 



bo(r) 



3^o(r)v 



(4) 



Combining Eqs. (||) and (|j), and using the fact that 
Vno(r) is proportional to Vt/o(r), we obtain for the 
equation of motion for v(r, i), 

|? = |^V[V.v] + V[v.f] + |[V.v]f. (5) 

at 3TO?io(r) 3 

Equation (^) has previously been derived by Griffin et al. 
H using kinetic theory. 

In our present discussion we assume that the potential 
is axially symmetric, 



U {r) = -mw%(x 2 +y 2 + \z 2 ), 



(6) 



where uq is the frequency of the trap in the x — y plane; 
A = lu^/loq - where to z is the frequency along the z-axis 
- expresses the anisotropy of the trap. Most experiments 
on oscillations have been done in such traps, but our 
results may be generalized to the case of traps with no 
axis of symmetry. As shown in Ref. 0, in the lowest 
modes, which in a spherical trap correspond to monopolc 
and quadrupole vibrations, the velocity field has the form 
v(r, t) = v(r)cos(wi), where v(r) = a(xic + yy) + bzz. 
Here a and b are constants, and the frequencies, to, of the 
two modes are 



1 



4A + 5±(16A 2 -32A + 25) 



1/2 



(7) 



and b/a = 3w 2 /2cj 2 — 5. 

We turn now to the damping of these modes. We adopt 
the standard approach of evaluating the rate of change of 
the mechanical energy, -Emech, associated with the mode, 
which is given by M 



E mcch = - I -\VT\ 2 dr 



C(V • v) 2 dr 



dvj 
dx k 



dvk 



dr, 



(8) 



where re is the thermal conductivity, r/ is the first, or 
shear, viscosity, £ is the second, or bulk, viscosity and 
T is the temperature. Because the system is inhomoge- 
neous, the transport coefficients are generally spatially 
dependent. Next we show that the shear viscosity is 
the only source of damping of the low-lying modes de- 
scribed above. The contribution from thermal conduc- 
tion vanishes because there are no temperature gradients 
for these modes, and the contribution from the bulk vis- 
cosity vanishes because £ vanishes. To demonstrate the 
absence of temperature gradients, we observe that for 
the modes under consideration, which have the velocity 
field given above, V • v = (2a + 6) cos(a;£) is indepen- 
dent of position. Thus, from the continuity equation, 
dn/dt + nV • v = 0, it follows that (dn/dt)/n is also 
constant. For an adiabatic process in a free monatomic 
gas, n oc T 3 / 2 , and therefore the deviation of the tem- 
perature from its equilibrium value is spatially indepen- 
dent. Consequently there are no temperature gradients 
generated by the mode, and hence no dissipation due to 
thermal conduction. As to the second viscosity, this van- 
ishes because a spatially-homogeneous, non-relativistic 
monatomic gas in equilibrium subjected to a slow uni- 
form dilation (r — » vr) remains in equilibrium, but at a 
different temperature. This result is independent of the 
statistics of the atoms and of the degree of degeneracy, 
and it is discussed in Ref . (To) . Inserting into Eq. (JsJ) the 
expression for the velocity field, we arrive at the follow- 
ing simple expression for the time-average of the rate of 
loss of mechanical energy: 



(-^rnech) 



bf / ?7(r)dr. 



(9) 



Equation (g|) implies that (-Emcch) vanishes for the 
monopole mode for the isotropic case (A = 1). 

The next task is to calculate the first viscosity. At the 
low energies of interest in experiments, the scattering of 
two atoms is purely s-wave, and the total cross section is 
a — 87ra 2 cat , where a sca t is the scattering length. We con- 
sider the viscosity in the classical limit, since we expect 
the classical limit to be quantitatively accurate, even at 
T very close to T c , from comparison to the known effects 
of degeneracy on the heat capacity. The viscosity has the 
general form pOl 



V = c v 



(mkT) 



1/2 



(10) 



A simple relaxation-time approximation, with a scatter- 
ing rate equal to no(r)crv, where v — (2e/m) 1 ' 2 is the 
particle velocity and e is the single-particle energy, leads 
to the result C r) = 2 7 / 2 /(15tt 1 / 2 ) « 0.426, while a varia- 
tional calculation gives C, ; = 57r 1 / 2 /2 4 w 0.554. 

An important feature of the expression for the viscosity 
is its independence of the particle density. Consequently, 
the integral in Eq. (||) formally diverges at large dis- 
tances from the center of the trap. The reason for this 
is that hydrodynamics fails, because in the outer parts 
of the cloud, particle mean free paths are too long for 
hydrodynamics to be applicable. To solve this problem, 
one should treat the outer parts of the cloud using kinetic 
theory, rather than hydrodynamics, but to obtain a first 
estimate of the effects we shall assume that the hydro- 
dynamic description holds out to a distance such that an 
atom incident from outside the cloud has a probability of 
no more than 1 /e of not suffering a collision with another 
atom. In the parlance of radiative transport, this corre- 
sponds to an optical depth of unity. Mathematically, this 
condition is 



1 



T 0,: 



ds 



(ii) 



where the local mean free path is ^(r) = [no^er]" 1 , and 
ro, s is the cutoff, which depends on direction. The in- 
tegral in Eq. ( |ll| ) is to be performed along the path for 
which the density gradient is steepest, that is along VE/o, 
and ds is the corresponding line clement. For a classical 
distribution, no(r) = n(0)e~ Uo ^ T ^ kBT , where n(0) is the 
density at the center. Since dUo — |V£/o|<is, Eq. jll| ) 
can be written as 



1 w n(0)ak B T- 



1 



vcr (r 0t ,)| 



-U (T ,s)/k B T 



(12) 



For large values of the dimensionless parameter 
X^ 2 Na 2 cat miu 2 /k B T, Eq. (pi) gives 



k B T 



1 



' 0,s 



moj 2 sin 2 9 + A cos 2 9 



lnr s (0), 



(13) 



where 9 is the angle of Yq :S with respect to the z-axis. 
The dimensionless quantity t s (9) is essentially the total 
optical depth at the center of the cloud 



m^r^' 2 



sin 2 9 + A cos 2 9 
sin 2 9 + A 2 cos 2 9 



1/2 



where 



t = crn(0) 



( k B T 
\ 2toAwq 



1/2 



(14) 



(15) 



The density at the center of the cloud is n(0) = 
N\ 1 / 2 (muj 2 /2Trk B T) 3 / 2 . The volume of the atomic cloud 



2 



is given by integrating 7q /3, Eq. (|l3|), over the solid an- 
gle. Equation (|9|) with the variational estimate for the 
viscosity (C v — 5tt 1 / 2 /2 4 ) can then be written as 



IE \~ 57:1/2 (a b\* {kBT)2 

The function /(A, r ) is defined as 
dx 



/(A, T ). (16) 



/(A, to) 



i [(A - l)x 2 + l]3/a 

(A-l)a; 2 + l 



In r A 



1/2 



(A 2 



lk 2 



1/2N 



3/2 



(17) 



We expect the leading term (oc [lnr ] 3 / 2 ) to be asymp- 
totically exact for large to , but the value of the cutoff in 
the logarithm will depend on the detailed kinetic- theory 
solution to the boundary-layer problem. 
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FIG. 1. We plot /(A, To) from Eq. @ as function of 
the number of particles N, for A = 8, a scat = 53 A, 
vo — 129 Hz, and T = 1.3 T c , which are the parameters 
in the experiment of Ref. [6]. Asymptotically, for m > 1, 
/(A = 8, to) « 2A- 1 / 2 [ln(r A 1 / 2 )] 3 '' 2 - 0.53 [1h(t AV 2 )]V2. 



To calculate the damping rate we need to evaluate the 
mechanical energy of the cloud, S mec h- In a- n y oscilla- 
tor the time-average of -EWch is equal to the maximum 
kinetic energy in the mode, so 



(3 



1 



mcch/ 



mno(r)w 2 (r)(ir, 



and therefore, 



(E me ch) — N 



k R T 



2A 



(18) 



(19) 



We can now introduce the amplitude damping rate 
(as opposed to the energy damping rate, which is 



damp 



twice the amplitude damping rate), given by the absolute 
value of the ratio (£ ; mech)/2(£ ; mcch), 



57T 1 / 2 k B T 1 (a-b) 2 
2 9 / 2 3 2 mwoaLt N (a 2 + 6 2 /2A) 



/(A, to). (20) 



In order to elucidate the origin of the damping of the 
modes given by Eq. (|7|), we remark that when the vis- 
cosity is taken to be constant, the non-equilibrium stress 
tensor is constant in space, and therefore does not give 
rise to any force on a fluid clement. The damping is 
caused by the entropy generated by the non-equilibrium 
energy current density, which has a constant, non-zero 
divergence. 

Let us now turn to the experiments on sound propa- 
gation that have been performed to date. A necessary 
condition for hydrodynamics to be applicable is that the 
mean free path be small compared with the characteris- 
tic length scale of the mode. For the low modes we are 
studying here, this is equivalent to the requirement that 
the optical depth t s (9) be large compared with unity in 
all directions. For the MIT experiment |Q, for which 
A < 1, this condition implies that A 1 / 2 ^ 3> 1. In terms 
of the characteristic lengths, R z — (2k B T /muj 2 ) 1 / 2 and 
R± = (2/csT/mwo) 1 ^ 2 ! which measure the spatial extent 
of the cloud along the z-axis and perpendicular to it, the 
condition A 1 / 2 ^ ^> 1 is seen to be equivalent to the re- 
quirement R± 3> 1(0), where Z(0) is the mean free path 
in the center of the cloud. For the parameters of this 
experiment Q| , the minimum number of particles that is 
required for this condition to be satisfied is w 3 x 10 9 , 
whereas the number of particles at T = 2 T c is found to 
be N ~ 5 x 10 7 pJ| . For the parameters in the JILA 
experiment the minimum value of t s (6) is To, since 
A > 1. In this case the condition t ^> 1 is equivalent 
to R z 3> ^(0). Therefore, for the JILA experiment, the 
minimum number of particles required to attain hydro- 
dynamic conditions is « 1 x 10 7 at T = 1.3 T c , which is 
to be compared with the experimental value N m 8 x 10 4 
at T w 1.3 T c [||. Figure 1 shows /(A, To) as function 
of N, with all the other parameters equal to the ones of 
Ref. S and T — 1.3 T c . The lowest value of N is chosen 
to be N = 1.1 x 10 7 , i.e., the one for which To = 1. 

Now we estimate for the MIT and JILA experiments 
the magnitude of the characteristic lengths R z and R±, 
and compare them with the mean free path in the cen- 
ter of the cloud. The transition temperature is obtained 
from k B T c = 0.94 iVV^^gw-) 1 ^. For the MIT experi- 
ment TV « 2.5 x 10 7 at T c 0, while uj w 2tt x 250 Hz 
and A m 5.8 x 10~ 3 [§. At T = 2 T c the particle number 
is 5 x 10 7 , resulting in R z i=s 380 /j,m and R± « 29 /im. 
The corresponding value of the central density n(0) = 
N/(n 3 ^R z R 2 L ) is w 2.9 x 10 13 cm" 3 . With the scattering 
cross section am 1.9 x 10~ 12 cm 2 (a scat « 28 A) we 
obtain 1(0) = ^(O)^]" 1 w 180 /im, which is less than R z , 
but much larger than R± . The other condition for hydro- 
dynamic behaviour is that the frequency of the mode be 
small compared with an average particle scattering rate. 
A single particle undergoes collisions at a rate no(r)cn'th ) 
where Vth — (SksT/trm) 1 / 2 is the average particle ve- 
locity. The average scattering rate in the cloud is thus 
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T scat = / nl(r)av th dr/ f n (r) dr = n(0)<r(k B T /ttto) 1 / 2 , 
since / ng(r) dr/ J n (r) dr = n(0)/2 3 / 2 . For the MIT ex- 
periment we therefore estimate wr sca t ~ 2.2 at T = 2T C . 

For the JILA experiment loo ks 2tt x 129 Hz, A = 8 
and the particle number at T c is N w 4 x 10 4 0. At 
T = 1.3 T c , where iV « 8 x 10 4 , we obtain R z « 4 ^m 
and 10 /im, while the central density is n(0) « 

3.7 x 10 13 cm -3 . The mean free path 1(0) is estimated 
to be 38 /im corresponding to a scattering cross section 
a w 7 x 10~ 12 cm 2 . Finally the dimensionless parameter 



WT SC at is ~ 19. We should mention that using Eq. (20) we 
find that for N = 1.1 x 10 7 , the amplitude damping times 
for the JILA experiment are « 8 ms and 47 ms at T = 
1.3 T c for the modes which correspond to monopole and 
quadrupole vibrations in a spherical trap, respectively. 

We thus conclude that even though in the MIT ex- 
periment 1(0) ^R z , there are as yet no experiments with 
which we can directly compare our calculation of the 
damping. To obtain a semi-quantitative description, we 
adopt a phenomenological interpolation formula for the 
frequency and the damping rate of the modes. This has 
the usual form, 



(21) 



characteristic of relaxation processes. Here u>c is the fre- 
quency of the mode in the colhsionless regime and loh i n 
the hydrodynamic regime, and r is a characteristic relax- 
ation time, which we anticipate will be of order the scat- 
tering time. Equation ( pl| ) gives the qualitatively correct 
limiting behaviour for lot 3> 1 and lot <C 1. The imag- 
inary part of the above equation gives for the damping 
rate 
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damp 



IK 

2 1 + (lot) 2 



(22) 



The amplitude damping time in the MIT experiment 
H is measured to be 80 ms at T — 2 T c . The frequency 
u c = 2v z « 38 Hz; Eq. (@) gives for v H w (12/5) 1 / 2 I / 3 « 
30 Hz. Using these numbers and Eq. (p2[), we can solve 
for wt to get the two solutions 3.6 and 3.6 -1 « 0.28. We 
regard the larger solution as being the physically relevant 
one, since it is close to our estimate of wr sca t « 2.2. The 
real part of Eq. (pl|), with lot — 3.6, implies that the 
oscillation frequency is w 37.5 Hz, which is consistent 
with the experimental value of 35 ± 4 Hz. For the JILA 
experiment Q, the amplitude damping time is ~ 50 ms 
at T w 1.3 T c . Since v c « 258 Hz and i/jj w 221 Hz, 
Eq. ( p2] ) implies that wt w 10.6 or 10. 6" 1 . Again we 
consider the larger value to be the physical one, since it 
is close to our theoretical estimate above, wr scat w 19. 
We thus conclude that damping of the modes in clouds 
of bosons above T c is in good agreement with theoretical 
expectations. 
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